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ABSTRACT 

We analyze a sample of 30,000 nearby obscured AGNs with optical spectra from the Sloan 
Digital Sky Survey and mid-IR photometry from the Wide-field Infrared Survey Explorer. 
Our aim is to investigate the host galaxy properties in AGNs where the nuclear activity level 
is indicated by the mid-IR luminosity, and to compare the results to previous studies, which 
have used the [Olll] emission line as the main black hole activity indicator. We first carry out 
a systematic study of how the mid-IR colours of AGN hosts vary as a function of age-sensitive 
4000 A break strength and optical "Eddington parameter" (L[Oiii]/M B h)- We find that the 
[3.4] - [4.6] micron colour has weak dependence on D n (4000), but strong dependence on 
-L[Oiii]/Mbh- We then use a "pair-matching" technique to subtract the 4.6 micron stellar 
emission contributed by the host galaxy. Intrinsic 4.6 micron AGN luminosities can be recov- 
ered for most Seyferts, but only statistically for LINERs. By combining our sample of Seyferts 
with a sample of type 1 AGN and quasars at z < 0.7 from the SDSS, we show that the [Olll] 
and 4.6 micron luminosities correlate roughly linearly over 4 orders of magnitude, but with 
substantial scatter. We also compare the partition functions of the total integrated 4.6 micron 
and [Olll] line luminosities from Seyferts as a function of a variety of host galaxy properties, 
finding that they are identical. If we include the total 4.6 micron flux from the LINER popu- 
lation, we find that more of the total [Olll] luminosity arises from bulge-dominated galaxies 
with old stellar populations. This supports the hypothesis that an abundant supply of gas is a 
pre-requisite for the formation of an obscuring "torus". 

Key words: galaxies: active; galaxies, nuclei; infrared: galaxies. 



1 INTRODUCTION 

The energy output of active galactic nuclei (AGNs) is thought to be 
a good probe of black hole (BH) growth history in the Universe. In 
principle, observing the radiation from the innermost part of the ac- 
cretion system yields a direct estimate of the total mass of gas swal- 
lowed by the central black hole. Luminous type 1 quasars clearly 
exhibit a featureless continuum which originates from the hot ac- 
cretion disk. However, in many type 2 AGNs it is not possible to 
observe the accretion disk because the radiation is absorbed by the 
surrounding gas. These AGNs lack the power-law continuum and 
broad optical emission lines from the accretion disk. Instead, the 
spectral energy distribution (SED) is dominated by stellar emis- 
sion from the host galaxy. The presence of an actively accreting 
black hole is manifested by high ionization narrow emission lines 
arising from gas a few hundred parsecs away from the black hole 
that is bein g irradiated by the accretion disk. The u nification model 
of AGNs JAntonu cci 1993; Urrv & Padovanill 19951) postulates that 
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an anistrotropically distributed collection of absorbing clouds (of- 
ten referred to as the "torus") can lead to a natural explanation for 
varying types of AGNs with very different SEDs. 

The fact that the accretion disk is obscured in type 2 AGNs 
means that we can study the host galaxy due to greatly reduced 
contrast between core radiation and outer stellar emission. There 
have been longstanding efforts to answer key questions, such as 
what kind of galaxies host A GN and what tri g gers the accre- 
tion onto the black hole (e.g.. lHeckmanlll980l : So et alj|l997l : 
McLure et al. 1999). By investigating a large sample of type 2 
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AGNs, Kauffmann et al. 1 2003c) found that local AGNs are mostly 
hosted by galaxies with stellar masses greater than 10 10 Mq and 
with stellar surface densities greater than 3 x 10 8 M©. Black hole 
growth is currently occurring in low mass black holes located 
in low mass bulges, w hich are also still forming stars at present 
JHeckman et alj 12004). The ratio between the current black hole 
growth rate by accretion and bulge growth rate by star formation is 
~0.001, consistent with the observed Mbh — M^uiae relation (e.g., 
iMarconi & Huntll2003l ; lHaring & Ri xll2004l) . 

Further inspection of the relation between star formation in the 
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host galaxy and the accretion rate onto the black hole suggests that 
there are two modes of AGN accretion ( IKauffmann & Heckmanl 
2009). The majority of AGN hosts have little star formation and old 
stellar populations. These AGNs may be fed by stellar winds from 
evolved stars. Their inferred accretion rates are lo w on average, but 
the du t y cycle for this kind of activi ty is high (Kau ffmann et al.l 
12003d : IKauffmann & Heckmaij|2009l) . It has been suggested that 
a substantial fraction of such objects may not be "true" AGN, be- 
cause the low ionization lines could be produced by radiation f rom 
evolved stars dCid Fernandes et al.l201ltlYan & Blantonll2012h . 

A small fraction of AG N hosts show clear recent cen- 
tral star formation ( Kauffman n et al. I l2003cl : IWild et al.l 120071 : 
iKauffmann et"alll2007l ; lLit]|201Ch . These AGNs have higher ac- 
cretion rates on average and occur in galaxies where black holes 
and bulges are growing simultaneously. It has been suggested 
that a plentiful cold gas supply is the common source for optical 
AGN activity and central sta r formation (Kauffmann et al. 2007; 
IKauffmann & Heckmanll 20091) . 

Obscuration makes the estimation of accretion rate difficult. 
Because the accretion disk continuum is hidden, it is necessary 
to find indirect indicators to represent the accretion power. Usu- 
ally these indicators measure the amount of "reprocessed" radia- 
tion. In the optical, emission lines from the narrow line region are 
comm only used, especially the [Olll]5007 line (e.g. lHeckman et al.l 
2004). The narrow line emission arises from gas that extends over 
much larger scales than the accretion disk; the emission region is 
observed to extend over scales of a few hundred parsecs in type 2 
AGNs. The hard X-ray emission from the corona surrounding the 
accretion disk is another popular indicator of black hole accretion. 
At longer wavelengths, the optical/ultraviolet radiation absorbed by 
the torus is re-emitted in mid-infrared, and the mid-IR luminosity 
is expected to be correlated with accretion rate. Some mid-IR high- 
ionization emission lines are also used as AGN activity indicators 
JDasvra et al]|2008l ; JMelendez et alj|2008l) . At even longer wave- 
lengths, radio emission from the jet is another possible indicator of 
black hole accretion rate. 

Generally speaking, except for lower luminosity FR-I type ra- 
dio AGNs, these ind i cators are correlated with each other (e.g. , 
Heckmanetal. 2005; Melendez et al. 2008; LaMassaetal. 2010; 
Honig et al. 2010; Ichikawa et al. 2012). However, clear discrep- 
ancies have also been reported, leading to the worry that any AGN 
survey carried out at only one wavelength may bias our understand- 
ing. Lower-luminosity FR-I type radio AGNs are clearly distinct 
from optica l AGNs because they occur in more massive galaxies 
JBest et alj2005T) . 

The AGN samples f rom deep X-ray surveys are considered to 
be close to complete (see Alexander & Hickoxll2012l for a review), 



yet iHeckman et alj ( 120051) raise the possibility that many [Om] 
bright AGNs are missed from hard X-ray surveys. The analysis of 
the cosmic X-ray background suggests the existence of a class of 
AGNs with very high colum n density along the line of sight, or s o 
called Compton-thick AGNs JGilli et alj2007l ; lTreister et al.l2009h . 
It is commonly agreed that this kind of objects are missed even by 
the deepest X-ray survey but (at least) some of them are identi- 
fied in optical/IR observatio ns (Heckman et al. 2005; P anessa et al.l 
120061 ; iMelendez et"al]|2008l ; IGoulding et alj|201lh . Similarly, op- 
tical identi fication tech niques may miss some strongly accreting 
black holes. iNetzer et all d2006l) and lTrouille & Bargerl (l2010h find 
that the [Om]/Lx ratio decreases with increasing X-ray luminos- 
ity, which implies that some X-ray AGNs ar e probably not identi- 
fied in optical surveys. IGoulding et al.N2010l) also argues that opti- 
cal surveys are missing the accretion around low mass black holes. 



Both X-ray and optical techniques are clearly affected by ex- 
tinction/absorption. However, short-wavelength radiation that is ab- 
sorbed by dust, will be re-radiated at IR wavelengths, so "missing" 
objects should be recovered in IR-selected surveys. Previous work 
has claimed th at the number of detected AGNs is gre atly increased 
in IR surveys JFiore et al.ll2008l;lDonlev et al.ll200*8h . Nevertheless 

it is also reported that IR-selected AGNs are biased against weak 
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and type 2 AGNs (Cardamone et al. 2008). Some X-ray selected 

AGNs have been found to have IR colours consistent with pure star- 
forming galaxies, suggesting that IR colour-selection techniques 
will miss accreting black ho les in host galaxies with strong star 
formation JBrusa et al. 2009). 

In order to understand fully how black holes grow, we need 
to study AGN at multiple wavelengths. Some efforts have been 
made to unify the view of AGNs from different wavelengths based 
on the data in some specific sky regio n with intensive multi- 
wavelength coverage. iHickox et al.l ( 120091) compare X-ray, IR and 
radio- selected AGNs using data from the AGN and Galaxy Evo- 
lution Survey (AGES). They find that the X-ray and IR-selected 
AGNs are similar in terms of host galaxy properties and cluster- 
ing strength, while radio AGNs are clearly more strongly clus- 
tered. Optically-faint AGNs have harder X-ray spectra than typ- 
ical unabsorbed AGNs. Differences between X-ray and IR AGN 
hosts are mainly attribu ted to the slight stellar/BH mass difference. 
iBongiorno et alj ( 120121) compile a catalogue of obscured and unob- 
scured AGNs from the Cosmic Evolution Survey (COSMOS) sky 
region. They perform detailed multi-wavelength SED decomposi- 
tion into stellar and AGN components, estimating the AGN accre- 
tion rate and host properties simultaneously. They confirm that the 
AGN fraction is higher in more massive galaxies, independent of 
accretion rate. 



These studies, however, are based on relatively small sam- 
ples. In this paper, we will perform an analysis of a large sam- 
ple of local obscured AGN s, based on the data from the Sloan 
Digital Sky Survey (SDSS, I York et all |2000|> and the W ide-field 
Infrared Survey Explorer ( WISE, Iwright et al.l 120101) . In sec- 
tion [2] we will describe the organization of the data from the 
SDSS/W75£ database and how we construct our local galaxy 
sample. In section [3] we will analyze the WISE colour proper- 
ties of local galaxies, including passive galaxies, star forming 
galaxies and AGNs. We explore how WISE colours change as 
a function of position in the 2-dimensional plane of 4000 A 
strength and black hole accretion rate, as quantified by the 
extinction-corrected [Om] line luminosity scaled by the black hole 
mass of the galaxy. This allows us to specify which colours are 
sensitive primarily to star formation, and which to AGN activity. 

We then compare the host galaxies of IR-selected AGN with 
those of optically-selected AGN in section [4] We do this by exam- 
ining how the total IR emission from the torus in nearby AGN is 
distributed across galaxies as a function of global properties such 
as stellar mass, stellar mass surface density, concentration, and spe- 
cific star formation rate (sSFR). We compare these IR emissiv- 
ity "distribution functions" to results obtained from integrating up 
the total [Om] emission from narrow-line regions (Heckman et al. 
2004). In section \5\ we examine AGNs that are clearly identified 
at IR wavelengths, but not at optical wavelengths, and ask whether 
and how they are different from the rest of the population. Finally 
we will discuss our results and give a brief summary in section[6] 
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Figure 1. The redshift versus stellar mass. The background grey-coded his- 
togram and the white contours show the distribution of the whole SDSS 
spectroscopic sample with a 14.5 < r < 17.6 cut on the r-band model 
magnitude. The coloured boxes indicate the different stellar mass bins and 
redshift ranges for the volume-limited subsamples we use. 
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Table 1. The stellar mass bins for the sub-samples that make up volume- 
limited sample. The corresponding redshift ranges and the number of 
sources in each bin are also listed. 



2 DATA AND SAMPLE 

2.1 Matching the SDSS spectroscopic and WISE photometric 
catalogues 

Our parent sample includes all galaxies from the MPA/JHU SDSS 
DR7 cataloguqj with magnitudes in the range 14.5 < r < 17. 6. 
The r > 14.5 cut is used to remove the nearby galaxies with large 
angular sizes, because both the SDSS and WISE pipeline photome- 
try will fail for highly extended sources. We limit our analysis to the 
redshift range 0.02 < z < 0.21, so we do not have to worry about 
evolutionary effects within our sample. Given the fact that the ma- 
jority of AGNs are hosted by massive galaxies (Kau ffmann et al.l 
l2003ch . we limit our study to galaxies with log(M*/M ) > 9.8. 

We first define a sample of galaxies that is "complete" in 
stellar mass. We do this by dividing galaxies with stellar masses 
log(A4»/M ) > 9.8 into bins of width 0.2 dex in log(M») and 
evaluating the redshift interval over which all such galaxies are de- 
tected in the SDSS spectroscopic sample. This is illustrated in detail 



1 http://www.mpa-garching.mpg.de/SDSS/ 



Any band 213789(98.85%) 21942(100.00%) 

W1+W2 209142(96.70%) 21330(97.21%) 

W4 56697 (26.22%) 21942 (100.00%) 

All bands 54324(25.12%) 21254(96.86%) 

Table 2. The detection rates in WISE bands. SI and S2 are samples we use 
in this paper (see text for detail). "W1+W2" means detection in both 3.4 
and 4.6 micron bands. "W4" means detection in 22 micron band. 



in Figure Q] The upper redshift limits show where the r-band flux 
limit means we can no longer detect all galaxies in the given stellar 
m ass range. Our cuts are simi lar to the sample definition adopted 
bv lvon der Linden et al.l ( 120 101) in their Figure 5. The lower redshift 
limits are imposed to avoid galaxies with angular sizes that are too 
large for the SDSS and WISE catalogue photometry to be reliable. 
There are a total of 216272 SDSS sources in these volume-limited 
samples. We list the sample details in Table[TJ 

The SDSS galaxies are matched to the WISE catalogue within 
a search radius of 3" from the optical position. Given the small as- 
trometry errors for both the SDSS and the WISE catalogues, the 
probability of mismatches is negligible. In order to get reliable 
flux measurements, we only use the fluxes with signal-to-noise ra- 
tio larger than 3. A large number of sources in our sample are 
extended sources. In this case, the default WISE pipeline profile- 
fitting photometry will underestimate the the total flux, so we use 
the total magnitudes derived from elliptical aperture photometry in- 
stead. The parameters of the elliptical apertures, such as axis ratios 
and position angles, are not derived from the WISE images them- 
selves, but are taken from the 2MASS Extended Source Catalog 
dSkrutskie et al.ll2006T) . For more detailed discussion of the ellip- 
tical aperture pho tometry, please see section IV.4.C and VI. 3. e of 
ICutrietal.lj2012l) . 

The WISE Vega magnitudes are converted to AB magnitudes 
and absolute fluxes. The SED shape affects the conversion from the 
total flux to m onochromatic flux, s o additional colour corrections 
are necessary. Iwright et al.l ( 120 101) tabulate the relation between 
spectral shape and colour correction factors. Here for simplicity, 
we only use the table entries for power law forms F oc v a , and 
we interpolate between the WISE colours provided in the table to 
derive our corrections, which are typically less than 3% for the 3.6, 
4.6 and 22 micron bands (the 12 micron band is very broad, so the 
corrections can be as high as 10%). 

In the whole sample, 213789 (98.9%) sources are detected by 
WISE in the 3.4 and 4.6 micron bands at the 3-<r level and above. 
Only 54324 (25.1%) sources have > 3<r detections in all WISE 
bands. The 22 micron band has the lowest detection rate. We thus 
contruct a subsample by adopting a 22 micron flux cut of 7 mjy. 
This flux level is where the WISE images with average number of 
sc an frames will rec over ~95% of "real" sources (see section VI. 5c 
of ICutri et alj|2012l for more details). There are 21942 sources in 
this subsample (hereafter we call the whole volume-limited sample 
as SI and this subsample as S2), and 21254 (96.9%) of them are 
detected in all WISE bands. We use S2 to study the colour distribu- 
tions of galaxies in section [3] In Table [2] we summarize the WISE 
detection rates of our various samples. 

We weight each galaxy by the inverse of V ma x, which is de- 
fined as the total volume within which the galaxy could be located 
and make it into our sample. Because SI is a sample selected by 
stellar mass, while S2 is a sample selected by stellar mass and 22 
micron flux, we use different weightings for S 1 and S2. Figure [2] 
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Figure 2. The normalized distributions of redshift, stellar mass, K-corrected 
g-i colour and r-band concentration index (defined as the ratio of the radii 
enclosing 90% and 50% of the total r-band light). The black and red lines 
are for samples SI and S2, respectively. 



shows the distribution of a number of different galaxy properties 
for galaxies in samples SI and S2. Comparing with SI, S2 galax- 
ies tend to have lower stellar masses, bluer g-i colours and lower 
concentrations. As we will show in the next section, the differences 
arise because passive galaxies are generally not detected at 22 mi- 
cron. We note that the shape of the stellar mass distribution of SI 
galaxies is consistent with the stellar mass function calculated by 
iBell et all ([2003) from 2MASS/SDSS data. 



2.2 Optical classification 

Here we use the class i cal [Nh]/Hq versus [OlH]/H/3 di- 
agnostics (Baldwin et al. 1981; Veille ux & Osterbrockl 1 19871; 
iKauffmann et al. 2003c; K ewlev et alj 2006) to classify the galax- 
ies as star-forming or AGN from their optical emission line 
measurements. The fluxes of the key emission lines (Ha, H/3, 
[OIII]A5007, [Nll]A6584 and [Sll]AA67 17,31) are directly taken 
from the MPA/JHU catalogue. However, the errors on the emission 
line measurements from the pipeline are usually underestimated. 
We follow the recommendations on the webpage of MPA/JHU cat- 
alogw0and multiply the errors by 2.473, 1.882, 1.566, 2.039 and 
1.621, respectively. To ensure the classifications are reliable, only 
the emission lines with S/N ^ 3 are considered. In SI the re are 
27755 sources above the line that IKauffmann et alj (2003cJ) sug- 
gest to separate AGNs from star-forming galaxies. We call them 

"optical AGNs" here after. 

Following Heck man et alj y004), we use the [Oni] line lu- 
minosity as an indicator of AGN activity /black hole accretion rate. 
We use the B aimer decrement to correct the [Qui] line for dust ex- 
tinction, adopting the reddening curve in lWild et al.l (120071) and an 
intrinsic Ha/H/3 ratio of 2.87 for star-forming galaxies and 3.1 for 
AGN (Osterbrock 1989). In order to estimate the [Oni] luminosity 
from the narrow-line region, we use the simple method suggested 
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All 


216272 (100.00%) 


21942(100.00%) 


AGN 


27755 (12.83%) 


8133 (37.07%) 


Strong AGN 


7613 (3.52%) 


4105(18.71%) 


Seyfert 


9776 (4.52%) 


4171 (19.01%) 


LINER 


16377 (7.57%) 


3873 (17.65%) 


SF 


23604 (10.91%) 


11077(50.48%) 


Passive 


49401 (22.84%) 


9 (0.04%) 


Non-AGN 


73005 (33.76%) 


1 1086 (50.52%) 


Ambiguous 


115512(53.41%) 


2723 (12.41%) 



Table 3. Sample SI and S2. The first column shows the names of optical 
classifications. The numbers and the fractions are listed. Please see the text 
for the details of the sample definition and the classification criterion. 



bv lKauffmann & Heckmanl (2009) to separate the total [Oni] lumi- 
nosity into AGN and star formation components. The fraction of 
AGN contribution is calculated based on the galaxy's position on 
the BPT diagram (see their Figure 3). 

In this work, "strong" AGNs are defined as optically- 
identified AGNs with [Qui] luminosities larger than 10 7 Lq. 
IKauffmann et alj (2fJ03cj) show that false BPT classification be- 
cause of dilution by emission from HII regions in the surrounding 
host galaxy falling within the SDSS fibre aperture is not important 
for AGNs with [Oni] luminosities greater than this valuqj. 

We also classify the AGNs into Seyfert galaxies and LIN- 
ERs (low ionization nuclear emission-line regions), according to 
the [S hI/Hq ratio using Function 7 in the paper by Kewlev et al.l 
(2006). A small fraction of AGNs are classified neither as Seyferts 
nor as LINERs, simply due to low signal-to-noise ratio of the [Sll] 
lines. 

We select a sample of non-AGN hosts ("non-AGN"), which 
are either star forming galaxies ("SF") or galaxies without de- 
tected emission lines (S/N < 2, "passive"). Passive galaxies are 
also required to have large concentration index R90/R50 > 2.6, 
high stellar mass surface density log(/i*/(MQ/kpc 2 )) > 8.5 and 
large 4000 A break strength D n (4000) > 1.6. These values are 
the points where sharp transitions from young st ar-forming galax- 
ies to o ld passive galaxies are observed to occur dKauffmann et alj 
l2003bl) . 

In table [3] we list the source numbers of each galaxy type. 
Because of the 22 micron luminosity cut, the fraction of passive 
galaxies in sample S2 is much smaller than in sample SI. We note 
alarge fraction of objects (115512, 53.4% of SI, "ambiguous") are 
not classified into any of the subclasses described above due to our 
strict criterion. We do not use them in our further analysis. 



2.3 SDSS quasars 

In our SI and S2 samples, only narrow-line AGNs are included. We 
have extracted a sample of low redshif t (z < 0.7) type 1 AGN from 
the SP SS DR7 quasar catalogue ( She n et alj|201 lUSchneider et alj 
120101) . The upper redshift limit is chosen to ensure that the [Om] 
line still falls in the SDSS spectrum. We match this sample to the 
WISE catalogue within a 3" matching radius. Our type 1 sample 
includes 3165 quasars, of which 3086 (97.5%) are detected in all 
WISE bands. Since the quasars are usually core-dominated, we use 



3 We carried out test similar to what Kauff mann et alj ^2003c) did by calcu- 
lating the fraction of strong AGN in narrow bins of stellar mass and redshift, 
and our results are similar to theirs. 
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Figure 3. The WISE monochromatic luminosity as a function of stellar mass 
for S2 non-AGN galaxies. The blue, green and red points are for different 
4000 A break strengths: 1.0 < £>„(4000) < 1.3, 1.3 < D„(4000) < 
1.5, 1.5 < -Dn(4000) < 1.8. Any data bin with source number lower 
than 20 is dropped. The error bar shows the error on the median value at a 
confidence level of 95%, estimated by bootstrapping within each data bin. 
The dotted lines are the best linear fits to the data, assuming a slope of 1 . 



WISE magnitudes based on PSF-fit photometry. The [Olll] line 
is corrected for extinction using the Balmer decrement measured 
from the narrow components of the Ha and H/3 lines. There are 
592 quasars with reliable narrow components measurements to do 
this extinction correction. We do K-correction to the WISE lumi- 
nosities and colours of the quasars, using the QSOl template from 
IPolletta et all fe007t) . 



3 MID-IR COLOURS OF LOCAL GALAXIES 

3.1 Stellar emission 

The WISE 3.4 and 4.6 micron bands are very similar to the Spitzer 
IRAC 3.6 and 4.5 micron bands. Previous studies show that in 
galaxies, the emission in this wavelength range is dominated by 
older sta rs, and luminosities are thus strongly correlated with stel - 
larmass dMeidt et al.ll2012l ; lEskew et al.ll2012l ; lHwang et alj2012h . 
We plot the L — M* relations of S2 inactive galaxies in Fig- 
ure [3] The four panels show luminosities evaluated at 3.4, 4.6, 
12 and 22 micron. In order to investigate the sensitivity of WISE 
luminosites to star formation, we split our sample into 3 differ- 
ent bins of 4000 A break strength: 1.0 < D„(4000) < 1.3; 
1.3 < £>„(4000) < 1.5; 1.5 < D„(4000) < 1.8. The 4000 A 
break strength, D„(4000), can be regarded as an indicator of the 
specific star formation rate of the galaxy averaged over a timescale 
of around a Gyr. Unlike emission line fluxes (such as Ha), the 4000 
A break is insensitive to extinction and to "contamination" from 
AGN. It allows a direct compari son of the stellar populations of 
AGN hosts and inactive galaxies JKauffmann et a l. 2003a). As can 
be seen from Figure|3] all the WISE luminosities have clear positive 
correlation with stellar mass, with slopes close to 1. With increas- 
ing wavelength, the difference between young and old galaxies be- 
comes larger. This confirms that the 3.4 and 4.6 micron bands are 



Figure 4. The distribution of SI non-AGN galaxies is displayed as grey 
scaled histogram background on the [3.4] - [4.6] colour versus 4000 A break 
plane. Single stellar population models are overploted as blue (sub-solar, Z 
= 0.008) and green (solar metallic! ties, Z = 0.02) curves. The circles and 
squares are from BC03 and CB07 respectively. The data points along the 
curve are from the templates with stellar age of 0.005, 0.025, 0.1, 0.29, 
0.64, 0.9, 1.4, 2.5, 5 and 1 1 Gyr. The model curves for each redshift bin are 
from the models convolved with WISE band filters at redshifts 0.03, 0.05, 
0.08, 0.12 and 0.17, respectively. 



indeed dominated by the light from old stars. Younger galaxies have 
lower mass-to-light ratios, probably due to additional contribution 
from thermally pulsing asymptotic giant branch (TP-AGB) stars. 
At 12 and 22 microns, the luminosity difference between galax- 
ies of the same stellar mass with different 4000 A break strengths 
becomes very large. In these bands the emission from old stars is 
no longer dominant. Instead, the dust emission makes a substantial 
contribution to the total flux. 

In Figure [4] we show the distribution of SI non-AGN galax- 
ies on the [3.4] - [4.6] colour versus 4000 A break (i.e. stellar 
age) plane. Since our aim in this section is to investigate the emis- 
sion from stars, we make use of the SI sample, which is not bi- 
ased against galaxies with no ongoing star formation. Rather than 
K-correcting the colours, we show results in 5 narrow redshift 
"slices". The non-AGN sources show a clear bimodal distribution 
on the colour-age plane, reflecting the star- forming and passive 
populations of nearby galaxies. On average, the young star- forming 
galaxies are redder than older galaxies. 

We compare our data to the stellar population synthesis mod- 
els of lBruzual& Chariot! J2003I) (BC03, circles) and an updated 
version of these model (CB07, squares). The major difference be- 
tween the two models is the treatment of TP-AGB stars. In CB07 
the dusty TP-AGB emission is more important and the models thus 
predict much redder IR colours at intermediate ages. In the plot, 
we show the predicted location of a "simple stellar population" 
(SSP) at a range of different ages, and for two different metallicities 
(see figure caption for details). The CB07 SSPs are clearly redder 
than the BC03 SSPs at stellar ages between ~200 Myr and ~2 
Gyr, forming a red "bump" on the low metallicity CB07 curve. The 
[3.4] - [4.6] colours of the galaxies with the lowest 4000 A break 
strengths are well matched to the lower-metallicity single stellar 
population colours from the CB07 models, but not from BC03. This 
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[12] - [22] 

Figure 5. The WISE colour-colour diagram. The distribution of S2 galax- 
ies is plotted as grey-scale background. The blue, green and red contours 
indicate the distributions of SF galaxies, weak AGNs and strong AGNs, 
respectively. The cyan contour is for the SDSS quasars. 



implies that the [3.4] - [4.6] colour difference between young and 
old galaxies is mainly due to emisson from TP-AGB stars. More 
detailed fits to models with more realistic star formation histories 
are required to draw more quantitative conclusions. We defer this 
to future work. 

None of the model curves predict [3.4] - [4.6] colours redder 
than ~ 0.7 and there appear to be very few normal galaxies with 
colours redder than this value in the data. In the next section, we 
analyze the mid-IR colours of AGN host galaxies. 



3.2 AGN host galaxies 

Figure [5] displays the distribution of S2 galaxies in the [3.4] - 
[4.6] versus [12] - [22] WISE colour-colour plane. As can be seen, 
there is a clear peak at [3.4] - [4.6] colours near zero and [12] - 
[22] colours around 2. We have indicated the locations of different 
galaxy sub-populations, as well as quasars, using coloured contours 
as described in the figure caption. Weak AGNs have very similar 
colour distribution as star-forming galaxies. The peak of the colour 
distribution of strong AGNs is shifted with respect to that of star- 
forming galaxies, but the overlap between the two populations is 
very substantial. Only the SDSS quasars have WISE colours that are 
clearly disjoint from those of normal galaxies. As the host galaxy 
contamination is small for these systems, quasar colours reflect a 
"pure" AGN SED profile. We thus conclude that the [3.4] - [4.6] 
colours of most nearby type 2 AGNs are strongly affected by emis- 
sion from stars. Only a minority of the strong type 2 AGNs have 
[3.4] - [4.6] and [12] - [22] colours that are similar to quasars. 

We now examine how the WISE colours of AGNs vary as a 
function of 4000 A break strength (i.e. amount of recent star forma- 
tion in the host galaxy) and as a function of our optically-defined 
accretion rate indicator based on the extinction-corre cted [Qui] line 
lumin osity. We adopt the Mbh — <J relation from iGultekin et al.l 
(2009) to calculate the BH mass: 
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Figure 6. The WISE colours as function of 4000 A break. The plots are for 
sample S2 galaxies with 0.07 < z < 0.11. Blue curves are for normal 
star-forming galaxies. The green, orange and red lines are for strong AGNs 
in different Eddington ratio bins: log(L[OlII]/M B H) < -0.56, -0.56 < 
log(L[Olll]/M BH ) < 0.01, log(L[Oni]/M BH ) > 0.01, respectively 
(all the values are in solar units). The error bar is the 1-cr error on the median 
value, catculated from bootstrap resampling each data bin. Any data bin 
with source number less than 20 is discarded. 



Throughout this paper, we will use the Eddington parameter 
L[Om]/AfBH as the optical indicator of the central BH accretion 
activity level. 

Figure [6] shows the four different WISE colours of strong 
AGNs as function of 4000 A break (i.e. host galaxy sSFR). Results 
are shown for AGN divided into 3 different ranges in Z/[Oiii]/Mbh 
(see figure caption). Results for star-forming galaxies are also plot- 
ted in blue for comparison. In order to minimize the scatter due to 
the fact that our colours have not been K-corrected, we only use the 
galaxies within a relatively small redshift range 0.07 < z < 0.11. 
The redshift cut we use includes galaxies in the stellar mass range 
10.4 < log(M*/Mo) < 11.6 (see Figure[T}. Our conclusions do 
not change if the analysis is done in other redshift and stellar mass 
ranges. 

WISE colours are redder at higher sSFR. The [3.4]-[12] WISE 
colour is most sensitive to star formation. There is no difference 
between AGNs and star-forming galaxies in the top right panel 
of Figure [6] indicating that the [3.4]-[12] colour is insensitive to 
AGN activity. In other WISE colours, the AGN contribution is more 
prominent and stronger AGNs are redder at any fixed D n (4000) . 
The [3.4]-[4.6] and [12]-[22] colours are most sensitive to accretion 
rate and least sensitive to star formation. This is consistent with the 
observed peak shift between strong AGNs, weak AGNs and star- 
forming galaxies shown in Figure [5] 

In the following sections, we drop the 12 and 22 micron bands 
from our analysis, and focus on the 4.6 micron luminosity as our 
main indicator of AGN activity at mid-IR wavelengths. Although 
we see in Figure [6] that the difference in colour between different 
L[Oih]/Mbh bins is in fact larger for [12] - [22] than for [3.4] - 
[4.6], the loss in sample size incurred by requiring our AGN sam- 
ples to be complete at 22 microns is too large. All the plots in the 
rest of the paper are based on sample SI. 
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Figure 7. The number of objects as a function of [OlII] luminosity (left 
panel) and 4.6 micron luminosity (after host correction, right panel), nor- 
malized by the total number of S 1 AGNs. The solid and dashed lines are for 
all optically selected AGNs and Seyfert galaxies, respectively. In left panel, 
the blue curve is for all the objects in the subsample and the red curve is 
for the objects with positive 4.6 micron luminosities. In right panel, the red 
curve is the same population as left panel. The grey curve is the AGNs with 
negative 4.6 micron luminosities, binned by the absolute values of their lu- 
minosities. The error bars are estimated by bootstrapping within the whole 
SI AGN sample. 



4 MID-IR PROPERTIES OF LOCAL AGNS 

In this section, we will use the 4.6 micron luminosity, corrected for 
the contribution from stars, as a way to parametrize the IR prop- 
erties of the AGN in our sample. We will then compute how the 
total IR emissivity from AGN is distributed among galaxies with 
different masses, structural parameters and stellar populations. Fi- 
nally we will compare these distribution functions with similar ones 
computed for the total [OlII] emissivity. 
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Figure 8. 4.6 micron AGN luminosity versus [OlII] luminosity. The black 
dots and upper limits show results of stacked SI Seyferts in [OlII] lu- 
minosity bins. The solid line is a bisector linear fit to SI Seyferts with 
L[OlII]> 10 7 Lq. The fitting results are shown on the bottom-right cor- 
ner. The dashed line is a linear fit, assuming the 4.6 micron luminosity is 
proportional to [OlII] luminosity. The median value of the IR-to-[OlII] ratio 
is shown on top-left corner. The grey background contour shows the distri- 
bution of SI Seyferts only with positive 4.6 micron fluxes. The cyan dots 
are SDSS quasars. 



4.1 AGN IR luminosity 

In the previous section, we showed that the [3.4] - [4.6] colour is 
the best WISE indicator of AGN activity. This implies that the 4.6 
micron luminosity can be used as an estimate of the mid-IR lu- 
minosity that originates from the torus itself, provided that we are 
able to subtract the contribution that originates from stars. In this 
section, we devise a method for performing this subtraction. 

Because the SI sample is large, we find non-AGN control 
galaxies with properties that closely match those of the AGN host 
galaxies. We match the AGNs to non-AGN galaxies in stellar mass, 
4000 A break, redshift and stellar mass surface density. Matching 
in stellar mass and D n (4000) ensures that the control galaxies have 
the same total stellar masses and stellar population ages as the AGN 
hosts. Matching in both redshift and stellar surface mass density 
ensures that the 3" diameter SDSS fiber aperture samples the same 
fraction of the total light from the hosts in the two samples. In order 
to minimize the scatter, we use strict matching criteria: stellar mass 
within ±0.01 dex, 4000 A break within ±0.025, redshift within 
±0.02, and stellar mass surface density within ±0.2 dex. For every 
AGN we estimate the contribution of stars to the total 4.6 micron lu- 
minosity from the median 4.6 micron luminosity of all the matched 
control galaxies. Typically there are 8 non-AGN "neighbours" for 
each of AGN to give reasonable estimation of the non-AGN com- 
ponent. We then subtract this estimate of the non-AGN component 
from the observed 4.6 micron luminosity to get the "pure" AGN 
4.6 micron luminosity L4.6 Mm .AGN- This correction should be re- 
garded as a statistical one. 

In some cases, particularly when the AGN luminosity is low, 
the resulting flux will be negative. Figure [TJshows that a substantial 
fraction of objects have negative fluxes when L[Om]< 10 7 L© or 



L4.6fj.rn, agn < 3 x 10 8 Lq. We also estimate the uncertainty of 
individual object by calculating the scatter of the 4.6 micron lumi- 
nosities of its "neighbours" used for host subtraction. This yields 
typical uncertainty value of ~ 2 — 7 x 10 8 Lq. Below this lu- 
minosity level, the nuclear emission is poorly determined due to 
host contamination. Comparing with the whole AGN sample, the 
Seyfert galaxies are less affected due to their higher nuclear lumi- 
nosities. Seyferts with positive fluxes always account for a larger 
proportion even when summing up to the lowest luminosities bins. 
The 4.6 micron luminosities of Seyferts are much better recovered 
individually, allowing a direct comparison between IR and [OlII] 
luminosities. 

Figure [8] shows that there is a good correlation between the 
corrected 4.6 micron monochromatic "pure" AGN luminosity and 
the [OlII] luminosity for Seyferts, though the scatter is as large as 
~ 0.28 dex on average. We do not include the LINERs in the plot, 
because the 4.6 micron flux from the central source cannot be es- 
timated accurately due to host galaxy contamination. We stack the 
Seyferts in different [OlII] luminosity bins to reduce the uncertainty 
(black dots and upper limits). We limit the linear fitting to the ob- 
jects with L[0lll]> 3 x 10 6 L Q (66% of the Seyferts in volume- 
weighted number), because at lower luminosity level it is too noisy 
to recover nuclear emission even with the stacking technique. The 
linear fit to the bright narrow line Seyferts gives a correlation of 
L 4 .6 M m,AGN <x L[OIII]°- 96±0 ' 07 (solid line). If we assume the IR 
luminosity is proportional to [OlII] luminosity, then we get a me- 
dian IR-to-[OlIl] luminosity ratio of ~ 24 (dashed line). We note 
these values are only valid for bright objects. We also compare the 
results with SDSS quasars. The quasars extend the correlation a 
further 2 orders of magnitude in [OlII] and in IR luminosity. 
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Figure 9. The integrated [OlII] and IR emissivity distribution as a function of various AGN properties for Seyferts (see text for details). The blue histogram is 
for the [OlII] luminosity and the red histogram is for the 4.6 micron luminosity. 



4.2 Comparison of the distribution of total IR luminosity and 
total [Om] luminosity from AGN as a function of host 
galaxy properties 



iHeckman et alj (2004) investigated the integrated [OlII] luminos- 
ity from type 2 AGN binned up as a function of stellar mass, of 
stellar surface mass density, of concentration index and of 4000 A 
break strength Z? n (4000). They showed that most of the present- 
day accretion traced by the [OlII] line is taking place in galax- 
ies with young stellar ages (D n (4000) < 1.6), intermediate stel- 
lar masses (10 10 - few xlO 11 Mq), high surface mass densities 
(3 x 10 8 — 3 x 10 9 M /kpc 2 ), and intermediate concentrations 
(R90/R50 = 2.2 — 3.0). In this section we carry out the same 
exercise using the integrated 4.6 micron luminosity and compare 
the results with what is obtained for the integrated [OlII] luminos- 
ity. We note that we use 4.6 micron luminosities that are corrected 
for emission from stars and [OlII] luminosities that are corrected 



for extinction and for the contribution from star formation for this 
exercise. 

The results are shown in Figure|9] where blue histograms show 
the distribution of the total [OlII] emission from the Seyfert sample 
and the red histograms show the distribution of the total 4.6 micron 
luminosity as a function of a wide variety of different host galaxy 
parameters. From left to right, and from top to bottom, the host 
galaxy properties investigated in Figure |9]are the following: 

(i) The [OlII] line luminosity normalized by the black hole mass 
(Eddington parameter). 

(ii) The black hole mass estimated from the stellar velocity dis- 
persion. 

(iii) Stellar mass M*. 

(iv) Stellar mass surface density /i» . 

(v) The concentration index, defined as the ratio between 90% 
light radius and 50% light radius ratio R90/R50. 

(vi) Rest-frame g-i colour. 
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Figure 10. Similar to Figure[9]but emission from LINERs is now included. 
Results are shown only for two of the parameters: concentration index 
R90/R50 and 4000 A break D„(4000). 



(vii) 4000 A break, D„ (4000). 

(viii) Galaxy inclination estimated from the ratio of the major- 
to-minor axes, a/b. 

(ix) The Balmer decrement calculated from the ratio of Ha to 
H/3 line fluxes. 

Comparison of the red and blue histograms in Figure |9] in- 
dicates that the host galaxies of the Seyferts producing the inte- 
grated 4.6 micron luminosity and the integrated [Om] are identi- 
cal. This provides stro ng support for the standard Unified Model 
dUrrv & Padovanil 1 1995b . In Figure [To] we indicate the main 
changes to these results if we add in the [Olll] and 4.6 micron emis- 
sion from LINERS in the S 1 sample. The main shifts occur in the 
luminosity distributions as a function of concentration index (left) 
and 4000 A break strenth (right). The [Om] LINER emission be- 
comes skewed towards bulge-dominated systems with higher 4000 
A break strengths. This appears to support a picture where LIN- 
ERs cannot be integrated within the standard Unified Model either 
because the ionizing radiation arises from evolved stars or because 
a torus is absent. We again note, however, that in our SDSS/WISE 
sample, detection of low luminosity mid-infrared emission from the 
central source is impossible because of host galaxy contamination. 



5 INVESTIGATION OF THE SUB-POPULATION OF 
AGNS SELECTED USING MID-IR COLOURS 

In this section, we investigate the sub-population of galaxies iden- 
tified as AGN from their WISE photometry, but not from optical 
emission line diagnostics. We begin with a brief review of IR-only 
AGN selection methods. 



5.1 WISE IR AGN selection 

The original IR colour-colou r selection techn i ques were based 
on the Spitzer IRAC colours dLacy et alj|2004l : IStern et alj|2005l : 
iRichards et"alll2006| ; iDonlev et alj|2008h . The selection was then 
extended to WISE bands by synthesizing data from multiwave- 
length observatio ns and temp l ate stu dies I IAssef et al.ll2010f) . Later, 
using real data, [Stern et al.l (120121) adopted a simpler criterion: 
[3.4] — [4.6] > 0.8. We note that their selection criterion was tuned 
to AGN searches at higher redshifts. If we adopt this cut, we only 
find 435 sources in sample SI. 

Another method is IR p ower-law selection, which w as first 
applied to Spitzer IRAC data jAlonso-Herrero et al.ll2006l) and is 
considere d to be a stricter w ay to select AGN than colour-colour 
selection JDonlev et al]|2008l) . It is based on the fact that the AGN 




reduced x threshold 

Figure 11. The fraction of optically-identified AGNs in the population that 
meet our mid-IR power-law selection criterion a < —0.5 as a function 
of reduced \ 2 threshold. The solid circles and empty circles indicate strong 
and all optical AGNs, respectively. The error bars are simple Poisson count- 



g errors. 

IR selection 


ALL 


optical AGN 


optical strong AGN 


col 

pow 

col & pow 

all 


435 

503 
284 
654 


379(87.1%) 
434 (86.3%) 
255 (89.8%) 
558 (85.3%) 


358 (82.3%) 
394 (78.3%) 
243 (85.6%) 
509 (77.8%) 



Table 4. The number census of the IR selected AGN samples, "col" and 
"pow" are the [3.4] - [4.6] colour selection and the power law selection, 
respectively. Symbol "&" means the sources are selected by both methods, 
"all" means the sources selected by any of the IR selections. 

continuum in quasars is known to have a power-law form. In prac- 
tice, we make use of all four WISE bands and fit the broadband 
SED as follows: 



0.4 x Mab = —a x log(A/l (im) + c. 



(2) 



Here Mab is the monochromatic AB magnitude and A is the 
effective wavelength in each WISE bands. The free parameter a 
is spectral slope (/„ oc v a ). AGNs are required to have spectral 
slope that are sufficiently red (a < —0.5). The quality of the fit, 
i.e. the similarity of the SED shape to a pure power law, can be 
quantified by x 2 • A sample selected with looser \ 2 threshold will 
include more galaxies, but also be contaminated by more star form- 
ing galaxies. In FigureQT| we plot the fraction of optically-selected 
AGNs as a function of threshold in reduced x 2 statistics. Results 
are shown for all AGN (open symbols) and for strong AGN (filled 
symbols). As can be seen the fraction of optically-identified sources 
drops sharply above \ 2 ~ 1-5, particularly for strong AGN. We 
therefore select this as a threshold, which yields a sample of 503 
IR-selected power-law AGNs. 

Table [4] shows the number of the IR AGNs with different se- 
lection methods. In total, we find 654 IR AGN. The fraction of op- 
tically identified AGN is high (85.3%). We note that this is a much 
smaller number than could be identified optically. It is clear that IR 
AGN selection methods based on WISE colours will miss a large 
fraction of type 2 AGNs at low redshift. 
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Figure 12. Similar to Figure \5\ but the grey-scale background shows the 
location of strong optical AGNs. The green and red crosses indicate [3.4] - 
[4.6] colour and power-law selected AGNs. The red dashed line is the track 
of a pure power-law spectrum. 




Figure 13. Rest frame median SEDs of different subsamples of S2 galax- 
ies. All SEDs are normalized at ~1 micron (J band). Grey, green and red 
colours are for the whole S 1 sample, IR colour selected AGN and power 
law selected AGN, respectively. The error bar is the 1 a scatter within the 
bin. The four panels show SEDs for galaxy subsets with different optical 
classifications. 



5.2 The SEDs of AGN selected at mid-IR wavelengths 

In Figure[T2] we show the distributions of the IR AGN samples on 
the WISE colour-colour diagram as in Figure [3] Green crosses in- 
dicate sources selected by the simple [3.4] - [4.6] colour cut, while 
red crosses indicate power-law sources. The colours of IR AGNs 
are consistent with typical quasar SEDs and clearly avoid the locus 
of star-forming galaxies. 

Our photometric data covers 5 SDSS bands, 3 near-IR bands 
(2MASS bands, extracted from NYU Value-Added Galaxy Cata- 
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Figure 14. The number distributions of IR AGNs (red), strong opti- 
cal AGNs (green) and weak optical AGNs (blue). The properties shown 
here are [OlII] luminosity, Eddington parameter, 4000 A break and 
D( BPT). D(BPT) is the dis tance to the AGN/SF separation line defined 
by iKauffmann et al J (2003d ) on the BPT diagram. "Pure" AGNs have the 
largest D values. 



log produced bv lBlanton et al] (12005)) and 4 WISE bands. We use 
these bands to build a SED for each source in SI. We calculate the 
AB magnitudes in the rest frame, interpolated from neighbouring 
data points. The detection rate of WISE 12 and 22 micron bands is 
relatively low in sample SI, so here we use sample S2 instead. We 
note that if we use S 1 we get similar result. 

Figure[T3]shows that both IR selection methods lead to similar 
SED shapes. IR AGNs are similar to field galaxies in the optical and 
the near-IR, but clearly different beyond 3 micron. All galaxies, ex- 
cept the passive ones which have little dust emission, show a clear 
turnover at ~5 micron. This is the point where the dust emission 
starts to dominate the total radiation output for galaxies with ongo- 
ing star formation or nuclear activity. For IR AGNs, the turnover 
moves to shorter wavelength (~2-3 micron) and the SED is much 
flatter. This flattening strongly affects the [3.4] - [4.6] colours. We 
note there are 9 passive galaxies in S2, shown in top-right panel. 
The origin of their 22 micron fluxes is still unknown. We defer this 
to future work. 



5.3 Optical properties of AGN selected at mid-IR 
wavelengths 

Since both IR AGN selection methods lead to consistent SED 
shapes, we simply combine both IR AGN samples. In Figure [14] 
we compare IR-selected AGN with optically-identified weak and 
strong AGN. The properties we investigate are [OlII] luminosity, 
Eddington parameter, 4000 A break an d D(BPT), the dis t ance to 
the AGN/SF separation line defined by IKauffmann et al] J2003cl) 
on BPT diagram ("pure" AGNs have the largest D values). Unlike 
Figure[9] this plot shows fraction by number rather than fraction of 
the integrated IR or [OlII] emissivity. The main thing we learn is 
that AGNs selected at low redshifts using mid-IR colours have sim- 
ilar [OlII] luminosities and 4000 A break strengths as strong optical 
AGN, but have values of the Eddington parameters and the D pa- 
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rameter that are slightly higher. This bias arises because the mid-IR 
colours are much more sensitive to star formation in the host galaxy 
than the [Om]/H/3 and [Nll]/Ha ratios. As a result, only strongly 
accreting black holes with high Eddington parameters are selected 
by IR colour-based techniques. 



5.4 Optically-unidentified IR AGN 

Though the IR AGNs are found to be similar to strong optical 
AGNs, there is a small fraction (< 15%) of objects not identified 
as AGNs from their optical emission lines. This holds for both IR 
selection techniques. In total there are 96 (14.7%) IR AGNs that 
are not classified as AGN in the optical. We call them "IR-only" 
AGNs for short. 

There are two reasons why optical identification might fail. 
One is because they are mis-classified as star-forming systems. The 
optical classification fails because of emission-line contamination 
from the host galaxy. The second possibility is that at least one of 
the four emission lines required for reliable BPT classification is 
not detected. 

We find 20% of IR-only AGNs are optically identified as star- 
forming galaxies (hereafter we call them SF-IR-only AGNs). Most 
(15 out of 19 sources) are relatively metal rich (log([Nll]/Ha) > 
—0.6) and are located at the border between the AGN and star- 
forming populations. There are also 4-metal poor SF-IR-only 
AGNs (log([Nll]/H«) < -0.6) that fall on the left side of the 
BPT diagram. Metal-poor AGNs are rare and occur in less mas- 
sive galaxies JGroves et al.ll2006l) . Three of these objects do not 
have red [3.4] - [4.6] or [12] - [22] colours, so it is difficult to 
judge whether these galaxies are true AGNs or not. Interestingly, 
one object has [3.4] - [4.6] = 1.1 and [12] - [22] = 7.2. It is 
clearly an AGN-dominated object with relatively low stellar mass 
log(M»/M0) = 9.83. 

The majority of our IR-only AGNs (77 objects, 80% of the 
population), are optically unidentified, because one or more emis- 
sion lines are not detected with sufficient signal-to-noise ratio. In 
most cases, it is the H/3 line measurement that has low S/N. We 
can estimate a lower limit on their [OlIl]/H/3 ratio that places them 
well into the Seyfert regime. In conclusion, as far as we can tell, 
IR-selected AGN without optical emission line classification do not 
constitute a special class of object. 



6 DISCUSSION 

In this paper, we have matched a large sample of SDSS galaxies 
with redshifts in the range 0.02 < z < 0.21 with mid-IR pho- 
tometry from the WISE survey. Our aim was to investigate the host 
galaxy properties of AGN by using the mid-IR luminosity as our 
AGN activity indicator, and then to compare the results to previ- 
ous studies, which have used the [Olll] line luminosity as the main 
diagnostic of AGN activity. 

As an AGN activity indicator, the [Olll] line luminosity has 
the advantage that it is relatively insensitive to contamination by 
ionized gas excited by young stars in the host galaxy (the [Olll] 
luminosity produced in high metallicity HII regions is known to be 
weak). However, because the bulk of the [Olll] emission excited 
by radiation from the accretion disk arises from gas located at dis- 
tances of a few hundred parsecs from the galaxy center, the [Olll] 
luminosity is a rather indirect indicator of current accretion onto 
the central black hole. In contrast, recent high-resolution observa- 
tions indicate that the scale of the IR-emitting "torus" around the 



black h ole might be no more than a few parsecs in radius (e.g., 
i Bock et alj|2000l: ISwain et al.ll2003l: fprieto et alj|2004l; Ijaffe et al.l 



|2004| ; iPackham et alj|2005i JElitzud [20051 : iMason et alj|2006h . It 
thus provides a probe of accretion onto the black hole on much 
smaller scales. However, a significant fraction of the total mid-IR 
emission from galaxies arises from stars. In the 3.4 and 4.6 mi- 
cron bands, emission from stars with ages greater than ~ 1 Gyr 
dominates, and at longer wavelengths emission from the dusty in- 
terstellar medium becomes important. 

In this paper, we use the 4000 A break strength, D n (4000), 
as our main probe of present-to-past averaged star formation in 
galaxies. We first carry out a systematic study of how the mid- 
IR colours of AGN hosts vary as a function of both D n (4000) 
and optical "Eddington parameter" (L[Oiii]/A/bh), finding that 
the [3.4] - [4.6] micron colour to have the weakest dependence 
on 7J n (4000), but strong dependence on L[Oiii]/A/bh- W e use a 
"pair-matching" technique introduced by Kauffmann et al. (2006) 
to statistically subtract the 4.6 micron stellar emission contributed 
by the host galaxies of the the AGN in our sample, by extracting 
samples of non-AGN with similar redshifts, stellar masses, sizes 
and 4000 A break strengths as the AGN host galaxies. We use these 
corrected 4.6 micron luminosities to parametrize the strength of the 
central torus emission for the AGN in our sample. We show that 
intrinsic 4.6 micron AGN luminosities can be recovered for most 
Seyferts, but only statistically for LINERs. 

By combining our sample of Seyferts with a sample of type 1 
AGN and quasars at z < 0.7 from the SDSS, we show that [Olll] 
and 4.6 micron luminosities correlate roughly linearly over 4 orders 
of magnitude in luminosity, except the low luminosity end. How- 
ever, there is substantial scatter in this relation. To gain further in- 
sight, we carry out a systematic comparison of how the host galaxy 
properties of AGN change if the nuclear luminosity is parametrized 
by 4.6 micron luminosity instead of [Olll] luminosity. We quantify 
this change using the partition function of the total integrated 4.6 
micron/[Oni] line luminosity from type 2 AGN as a function of a 
variety of host galaxy properties including stellar mass, structural 
properties such as stellar surface mass density and bulge-to-disk ra- 
tio, and indicators of stellar population age and ISM dust content. 

We find identical distributions of total 4.6 micron and [Olll] 
line luminosity for Seyfert gala xies, in strong supp ort of the 
standard Unified Model (see also lLaMassa et alj|2012h . If emis- 
sion from LINERs is included, the integrated [Olll] emission 
shifts towards bulge-dominated galaxies with older stellar pop- 
ulations. This implies that LINERs may not have torus or that 
their [Qui] emission is not connected with B H accretion (e.g. 
ICid Fernandes et al.ll201 It lYan & Blantonl l2012). or some combi- 
nation of both effects. We also note that if we divide our sample 
by optical Eddington parameter or 4.6 micron luminosity scaled by 
black hole mass and if we repeat the comparisons using the 25% 
of the emission coming from the IR and optical sources with the 
highest accretion rates, host galaxy properties are also identical. 

Finally, we note that we searched the entire SDSS spectro- 
scopic catalogue for AGN that could only be identified as such us- 
ing WISE photometry. We found a total of 96 such systems. A de- 
tailed analysis revealed that there was nothing special about these 
objects: in most of them, the S/N in the H/3 line was simply too 
low to allow a reliable BPT classification. One might be tempted 
to conclude, therefore, that no differences exist between the opti- 
cal and IR "view" of low redshift AGN. In a companion paper, we 
present new results on close pair counts and asymmtery in the stel- 
lar light distribution of AGN as a function of both [Olll] and 4.6 
micron luminosity. 
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